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Abstract. Several problems of trigonometric approximation on a hexagon 
and a triangle are studied using the discrete Fourier transform and orthogo- 
nal polynomials of two variables. A discrete Fourier analysis on the regular 
hexagon is developed in detail, from which the analysis on the triangle is 
deduced. The results include cubature formulas and interpolation on these 
domains. In particular, a trigonometric Lagrange interpolation on a triangle is 
shown to satisfy an explicit compact formula, which is equivalent to the poly- 
nomial interpolation on a planer region bounded by Steiner's hypocycloid. The 
Lebesgue constant of the interpolation is shown to be in the order of (logn) 2 . 
Furthermore, a Gauss cubature is established on the hypocycloid. 



1. Introduction 

Approximation by trigonometric polynomials, or equivalently complex exponen- 
tials, is at the root of approximation theory. Many problems central to approxima- 
tion theory were first studied for trigonometric polynomials on the real line or on 
the unit circle, as illustrated in the classical treatise of Zygmund [26]. Much of the 
advance in the theory of trigonometric approximation is due to the periodicity of the 
function. In the case of one variable, one usually works with 2it periodic functions. 
A straightforward extension to several variables is the tensor product type, where 
one works with functions that are 27r periodic in each of their variables. There 
are, however, other definitions of periodicity in several variables, most notably the 
periodicity defined by a lattice, which is a discrete subgroup defined by AZ d , where 
A is a nonsingular matrix, and the periodic function satisifes f(x + Ak) — f(x), for 
all k S 1 d . With such a periodicity, one works with exponentials or trigonometric 
functions of the form e l2,ra ' a: , where a and x are in proper sets of K d , not necessarily 
the usual trigonometric polynomials. 

If CI is a bounded open set that tiles M. d with the lattice L — AZ d , meaning that 
fl + A7L d = M. d , then a theorem of Fuglede [B] states that the family of exponen- 
tials {e 2ma ' x : a S L- 1 }, where L 1 - = A~ tx 7L d is the dual lattice of L, forms an 
orthonormal basis in L 2 (Q). In particular, this shows that it is appropriate to use 
such exponentials for approximation. Naturally one can study Fourier series on f2, 
defined by these exponentials, which is more or less the usual multivariate Fourier 
series with a change of variables from 7L d to A7L d . Lattices appear prominently in 
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various problems in combinatorics and analysis (see, for example, j2[|5]), the Fourier 
series associated with them are studied in the context of information science and 
physics, such as sampling theory [JJ[TT] and multivariate signal processing [3]. As 
far as we are aware, they have not been studied much in approximation theory until 
recently. Part of the reason may lie in the complexity of the set Q. In the simplest 
non-tensor product case, it is a regular hexagon region on the plane, not one of 
those regions that people are used to in classical analysis. It turns out, however, 
that the results on the hexagon are instrumental for the study of trigonometric 
functions on the triangle [5TJ [55] . The Fourier series on the hexagon were studied 
from scratch in |21j . and two families of functions analogues to those of sine and 
cosine but defined using three directions were identified and studied. In |22j the two 
generalized trigonometric functions were derived as solutions of the eigenfunctions 
of the Laplace equation on a triangle, and the Fourier scries in terms of them were 
studied. 

The purpose of the present paper is to study discrete Fourier analysis on the 
hexagon and on the triangle; in particular, we will derive results on cubature and 
Lagrange interpolation on these domains. We shall follow two approaches. The first 
one relies on the discrete Fourier transforms associated with the lattices, which leads 
naturally to an interpolation formula, and the result applies to a general lattice. We 
then apply the result to the hexagonal tiling. The symmetry group of the hexagon 
lattice is the reflection group A2 ■ The generalized cosine and sine are invariant and 
anti-invariant projections, respectively, of the basic exponential function under A2, 
and they can be studied as functions defined on a fundamental domain, which is 
one of the equilateral triangles inside the hexagon. This allows us to pass from 
the hexagon to the triangle. The discrete Fourier transform can be regarded as a 
cubature formula. In the case of a triangle, the trigonometric functions for which 
such a cubature holds is a graded algebra due to a connection to polynomials of 
two variables, which shows that the cubature is a Gaussian type. 

The second approach uses orthogonal polynomials in two variables. It turns out 
that the generalized cosine and sine functions were studied earlier by Koornwindcr 
in [9], who also started with eigenfunctions of the Laplace operator and identified 
them as orthogonal polynomials on the region bounded by Steiners hypocycloid, 
which he called Chebyshev polynomials of the first and the second kind, respectively. 
It is well known that orthogonal polynomials are closely related to cubature for- 
mulas and polynomial interpolation (see, for example, [H [TJl [5D1 [2H Hi] ) . Roughly 
speaking, a cubature formula of high precision exists if the set of its nodes is an 
algebraic variety of a polynomial ideal generated by (quasi)orthogonal polynomials, 
and the size of the variety is equal to the codimension of the ideal. Furthermore, 
such a cubature formula arises from an interpolation polynomial. The rich struc- 
ture of the generalized Chebyshev polynomials allow us to study their common 
zeros and establish the cubature and the polynomial interpolation. One particular 
result shows that Gaussian cubature formula exists for one integral on the region 
bounded by the hypocycloid. It should be mentioned that, unlike the case of one 
variable, Gaussian cubature formulas rarely exist in several variables. 

The most concrete result, in our view, is the trigonometric interpolation on the 
triangle. For a standard triangle A := {(x,y) : x > 0, y > 0, x + y < 1}, the set of 
interpolation points are X n :— {(— } *j) : < i < j < n} in A. We found a compact 
formula for the fundamental interpolation functions in terms of the elementary 
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trigonometric functions, which ensures that the interpolation can be computed 
efficiently. Moreover, the uniform operator norm (the Lebesgue constant) of the 
interpolation is shown to grow in the order of (logn) 2 as n goes to infinity. This 
is in sharp contrast to the algebraic polynomial interpolation on the set X n , which 
also exists and can be given by simple formulas, but has an undesirable convergence 
behavior (see, for example, [17] ). In fact, it is well known that equally spaced points 
are not good for algebraic polynomial interpolation. For example, a classical result 
shows that the polynomial interpolation of f(x) — \x\ on equally spaced points of 
[— 1, 1] diverges at all points other than —1,0, 1. 

The paper is organized as follows. In the next section we state the basic results 
on lattices and discrete Fourier transform, and establish the connection to the in- 
terpolation formula. Section 3 contains results on the regular hexagon, obtained 
from applying the general result of the previous section. Section 4 studies the gen- 
eralized sine and cosine functions, and the discrete Fourier analysis on the triangle. 
In particular, it contains the results on the trigonometric interpolation. The gener- 
alized Chebyshev polynomials and their connection to the orthogonal polynomials 
of two variables are discussed in Section 5, and used to establish the existence of 
Gaussian cubature, and an alternative way to derive interpolation polynomial. 

2. Continuous and discrete Fourier analysis with lattice 

2.1. Lattice and Fourier series. A lattice in R d is a discrete subgroup that 
contains d linearly independent vectors, 

(2.1) L := {fciffli + k 2 a 2 H + k d a d : ki £ Z, i = 1, 2, . . . , d} , 

where oi, . . . , a d are linearly independent column vectors in M. d . Let A be the d x d 
matrix whose columns are ai, . . . , a d - Then A is called a generator matrix of the 
lattice L. We can write L as La and a short notation for La is A7L d \ that is, 

L A = AZ d = {Ak:ke Z d } . 

Any cZ-dimensional lattice has a dual lattice L 1 - given by 

L 1 - = {x e R d : x ■ y E Z for all y € L}, 

where x ■ y denote the usual Euclidean inner product of x and y. The generator 
matrix of the dual lattice L\ is A~ tr := (A tr )~ 1 , where A tr denote the matrix 
transpose of A ([!]). Throughout this paper whenever we need to treat x € R d as 
a vector, we treat it as a column vector, so that x ■ y = x tr y. 
A bounded set f2 of M. d is said to tile M. d with the lattice L if 

/J Xn(x + a) = 1, for almost all x € M d , 

where xsi denotes the characteristic function of fl. We write this as + L = M. d . 
Tiling and Fourier analysis are closely related; see, for example, 6,8 . We will need 
the following result due to Fuglede [B]. Let (•, -)n denote the inner product in the 
space L 2 {Vt) defined by 
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Theorem 2.1. [6] Let O in M. d be a bounded open domain and L be a lattice ofR d . 
Then f2 + L = M. d if and only if {e 2ma x a G L } is an orthonormal basis with 
respect to the inner product 



The orthonormal property is defined with respect to normalized Lebesgue mea- 
sure on f2. If L = La, then the measure of is equal to |detA|. Furthermore, 
since L\ = A~ tr Z d , we can write a = A~ tr k for a S L\ and fc S Z d , so that 
a ■ x = (A~ tr k) ■ x — k tr A x. Hence the orthogonality in the theorem is 

The set f2 is called a spectral (fundamental) set for the lattice L. If L = La we 
also write fl = Qa- 

A function / £ L 1 (flA) can be expanded into a Fourier series 

'« ~ E ■* - T*br /„ 

which is just a change of variables away from the usual multiple Fourier series. The 
Fourier transform / of a function / e L 1 (R d ) and its inverse are defined by 

/(£) = / /:•<•;<' JT ' " - '/•'• and /(a) = / /(£)e 2 ™« -df. 
One important tool for us is the Poisson summation formula given by 

(2.4) £ /(, + Ak) = £ f(A^k)e^ kt ' A ^, 

k^L d 1 V n fc6Z d 

which holds for / G L 1 (flA) under the usual assumption on the convergence of 
the series in both sides. This formula has been used by various authors (see, for 
example, [IJ [7J [9l [19]). An immediate consequence of the Poisson summation 
formula is the following sampling theorem (see, for example, [71 ITT|). 

Proposition 2.2. Let be a spectral set of the lattice AZ d . Assume that f is 
supported on fl and f £ L 2 (fl). Then 

(2.5) f(x) = Y, f{A- u k)m n {x - A~«k) 
in L 2 (Q), where 

In fact, since / is supported on Ha, /(£) = J2kez d /(£ + Ak). Integrating over 
f2, the sampling equation follows from (|2.4p and Fourier inversion. We note that 

Vn(A- tr j) =<y 0ji , foralljeZ d , 

by (|2.3p . so that \I/n can be considered as a cardinal interpolation function. 
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2.2. Discrete Fourier analysis and interpolation. A function / denned on R rf 
is called a periodic function with respect to the lattice A7L d if 

f(x + Ak) = f(x) for all k G Z d . 

For a given lattice L, its spectral set is evidently not unique. The orthogonality 
(|2.3[) is independent of the choice of J7. This can be seen directly from the fact that 
the function x i— » e 27Tlkt ' A lx is periodic with respect to the lattice AZ d . 

In order to carry out a discrete Fourier analysis with respect to a lattice, we shall 
fix f2 such that f2 contains in its interior and we further require in this subsection 
that the tiling f2 + AZ d = R d holds pointwise and without overlapping. In other 
words, we require 

(2.7) ^2 Xa{x + Ak) = 1, Vx G R d , and (O + Ak) n (fi + Aj) = 0, if fc ^ j. 

For example, we can take = [— i, for the standard lattice Z d . 

Definition 2.3. Lei A and B be two nonsingular matrices such that D,a o-nd Qb 
satisfy (|2.7|) . Assume that all entries of N := B tr A are integers. Denote 

A N := {fc G Z d : B~ tr k G fl A }, and A m , := {fc G Z d : A~ tr k G fl B }. 

Two points a;,?/ G R d are said to be congruent with respect to the lattice AZ d , 
if x — y G AZ d , and we write x = y mod A 

Theorem 2.4. Let A, £> and N be as in Definition \2.3\ Then 



(2 8 n 1 e 2 7 rifc*'iv- 1 j = I !. */ k = mod ^ tr > 

|det(JV)| £^ |0, otherwise. 

Proof. Let / € L 1 (R d ) be such that Poisson summation formula holds. Denote 

$ A f(x) := J2 f(x + Ak). 

kei d 

Then the Poisson summation formula (12.411 shows that 



$ A f{B-«j) = r^rr J2 fiA^k^'"- 1 ^ J G Z d . 

Since R d = f^s + BZ d , we can write A~ tr fc = xb + Bm for every fc G Z d with 
is £ fig and m G Z d . The fact that iV has integer entries shows that j4 tr iE_B = 
I G Z d . That is, every fc G Z d can be written as fc = / + N tr m with / G A^tr and 
m G Z d . It is easy to see that e 2, " fe N J = e 2 ™ 1 J under this decomposition. 
Consequently, 



|det(A)| 2- $ f[A l)e 

Similarly, we have 

(2.10) * B f(A-«l) = — * - £ * A f(B- tr 3)e 

1 v n jeAjv 
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Substituting (|2~9|) into (12.10|) leads to the identity 

1 1 feeA" iGA„ 



that holds for / in, say, the Schwartz class of functions, from which (|2.8|) follows 
immediately. □ 



Theorem 2.5. Let A,B and N be as in Definition \2.3\ Define 

jeA N 

/or /, g in C{SIa), the space of continuous junctions on D,a- Then 

(2.H) (f,9)n A = (f,g) N 

for all /, g in the finite dimensional subspace 

H N := span j^fc : 4> k {x) = e 2 ™^'^ 1 ^ fc e \ N A . 

Proof. By the definition of (-,-}n, equation (12. 8|) implies that {4>i,4>k)N = <5;,fc, 
l,k € N tr , which shows, by ([23]) . that (j2~TTj) holds for / = </>; and .g = fa- Since 
both inner products are sesquilinear, the stated result follows. □ 

If A is a subset of Z d , we denote by #A the number of points in A. Evidently 
the dimension oi7i N is #A N t,. 

Let T N f denote the Fourier expansion of / S C^a) m Hjv with respect to the 
inner product (•, •). Then 

1 v n k£A N jeAjytr 

Hence, analogous to the sampling theorem in Proposition ^. 21 we have 

(2.12) l N f(x) = ]T f(B- tr k)^ B (x - B- tr k), f e C(O0, 
feeAN 

where 

Theorem 2.6. Let A, B and N be as in Definition \2.3\ . Then T^f is the unique 
interpolation operator on N inTi^ ; that is, 

T N f(B-^j) = f(B~ tr j), Vj S A N . 

In particular, #Ajv = #Ajytr = | det(iV)|. Furthermore, the fundamental interpola- 
tion function satisfies 

(2-14) *£ B (x)=J2*n B (x + Aj). 
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Proof. Using (|2.8j) with N in place of N tr gives immediately that \I/q b (B tr fc) = 6k,a 
for k 6 Ajv, so that the interpolation holds. This also shows that {^q (2; — 



B tr fc) : fc S An} is linearly independent. The equation (|2.8p with fc = shows 
immediately that #A^v = |det(7V)|. Similarly, with N replaced by N tr , we have 
#Ajv = | det(7V tr )|. This shows that #A N = \ det(N)\ = #A N t,. Furthermore, 
let M = (<fik(B~ tr j))j t k£A N be the interpolation matrix, then the equation (|2.8p 
shows that M is a unitary matrix, which shows in particular that M is invertible. 
Consequently the interpolation by H n on points in Ajy is unique. 



Finally, since V&f 



| det (B) | 1 xn B , we have 



1 



j<£Z d 



\det(N)\ 



^^■(x)xo B (A- tr j) 



0j(a:)*n fl (A- tr j), 



from which (|2. 14)) follows from the Poisson summation formula (|2.4[) . 



□ 



The result in this section applies to fairly general lattices. In the following section 
we apply it to the hexagon lattice in two dimension. We hope to report results for 
higher dimensional lattices in a future work. 

3. Discrete Fourier analysis on the hexagon 

3.1. Hexagon lattice and Fourier analysis. The generator matrix and the spec- 
tral set of the hexagon lattice is given by 

\-l 2, 

The strict inequality in the definition of £Ih comes from our assumption in (|2.7[) . 



H 



]V = [(X1,X 2 ) : -1 < X 2 , ^Xl ± \x 2 < l} 



Figure 1. The lattice L H - 



As shown in [21] , it is more convenient to use homogeneous coordinates (t\ , t 2 , ts) 
that satisfies t\ + t 2 + t 3 = 0. We define 

^/3xi , , x 2 \[2>X\ 

Y 



(3.1) 



u=- x 4 



h = x 2 , t 3 = — — - 



8 



HUIYUAN LI, JIACHANG SUN, AND YUAN XU 



Under these homogeneous coordinates, the hexagon fin becomes 

fi = {(i1.t2.t3) eR 3 : -i<h,t 2 ,-t 3 <l;t 1 +t 2 + t 3 = 0}, 

which is the intersection of the plane t\ + t 2 +tj, = with the cube [—1, l] 3 , as shown 
in Figure 2. Later in the paper we shall depict the hexagon as a two dimensional 
figure, but still refer to its points by their homogeneous coordinates, as shown by 
the right hand side figure of Figure 2. 




FIGURE 2. The hexagon in homogeneous coordinates. 



For convenience, we adopt the convention of using bold letters, such as t, to 
denote points in the space 

R 3 H : = {t = (h,t 2 , tg) G R 3 : h + t 2 + t 3 = 0}. 

In other words, bold letters such as t stand for homogeneous coordinates. Further- 
more, we introduce the notation 

H := Z 3 n R% = {k = (fci, k 2 ,h) G 1? : k x + k 2 + fc 3 = 0} . 

The inner product on the hexagon under homogeneous coordinates becomes 

(/. 9)h = 777-7 / f(x 1 ,x 2 )g(x 1 ,x 2 )dx 1 dx 2 = rp-r I f(t)g(t)dt, 
Jn H Ju 

where |f2| denotes the area of Q. 

To apply the general results in the previous section, we choose A = H and 
B = ^H, where n is a positive integer. Then 



N = B tr A = —H T H 
2 



2n —n 
—n 2n 



has integer entries. Note that N is a symmetric matrix so that A^r = Ajyw. Using 
the fact that B~ tr j = ^H~ tr j = i + Ta^ 2 '-? 2 )' ^ * s n0 ^ nar d to see that 

B~ tr j G Q A , or equivalently j G A N , becomes j = (J1J2, -ji - h) G H„, where 

H„ := {j G H : -n < j x ,j 2 , -j 3 < n} and #H„ = in 2 . 
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In other words, An in the previous section becomes H„ in homogeneous coordinates. 
The fact that #H n = 3n 2 can be easily verified. Moreover, under the change of 
variables x = [x\, x 2 ) 1— > (ti, t 2 ,t 3 ) hi (|3.ip . it is easy to see that 

(3.2) (k 1 ,k 2 )H- 1 x = k 2 , -fei - k 2 )t = |k • t, 

so that e 27rlfctr - ff * = e^ k t . Therefore, introducing the notation 

0j(t) :=e^ ,rt , j€H, 
the orthogonality relation (|2.3p becomes 

(3.3) (<£k,<Aj)ff = <5k,j, kJeH. 

The finite dimensional space Hn in the previous section becomes 

TL n :— span{</)j : j G H„} , and dim7Y„ = 3n 2 . 

Under the homogeneous coordinates (|3.1[) . x = y (mod if) becomes t = s 
mod 3, where we define 

t = s mod 3 t\ — s\ = t 2 — s% = t-j — S3 mod 3. 

We call a function / H-periodic if /(t) = /(t + j) whenever j = mod 3. Since 
j, k £ H implies that 2j • k = (ji — j 2 )(ki — k 2 ) + 3j3&;3, we see that </>j is H-periodic. 
The following lemma is convenient in practice. 

Lemma 3.1. Let e\ = (2,-1,-1), e 2 = (-1,2,-1), s 3 = (-1,-1,2). Then a 
function f(t) is H-periodic if and only if 

(3.4) /(t + j £s: ) = /(t), jeZ, i = 1,2,3. 

3.2. Discrete Fourier analysis on the regular hexagon. Using the set-up in 
the previous subsection, Theorem 12.51 becomes, in homogeneous coordinates, the 
following: 

Proposition 3.2. For n > 0, define 

(3-5) (/.i4EM /,<?eC(Q). 

JSH„ 

Then 

{f,9)H = (f,g) n , f,geH n . 

The point set H„, hence the inner product (•,•)«., is not symmetric on f2 in the 
sense that it contains only part of the points on the boundary, as shown by the left 
hand figure in Figure 3. We denote by H* the symmetric point set on f2, 

H; := {j G H : -n < ji,j 2 ,j 3 < n} , and #H; = 3n 2 + 3n + 1. 

The set H* is shown by the right hand figure in Figure 3. 

Using periodicity, we can show that the inner product (•, •) is equivalent to a 
symmetric discrete inner product based on H* . Let 

H° :={j£M:-n<j 1 ,j2,h<n} 

denote the set of interior points in H* and H„. Define 

(3-6) (/,<?>; :=^E <f°/(£)5(J), f>9 e C(fi), 
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(n O n) ( n.O.nl 



Figure 3. The set HL and the set 



where 

!1, j £ H°, (inner points), 

i, jei;\i;,jiW3^0, (edge points), 

i, jeH* \H£, jd2h = 0, (vertices), 

in which we also include the positions of the points; the only one that needs an 
explanation is the "edge points" , which are the points located on the edges but not 
on the vertices of the hexagon. In analogy to Proposition 13. 21 we have the following 
theorem. 

Theorem 3.3. For n > 0, 

(f,a)H = (f,g) n = (f,g)* n , f,geH n . 

Proof. If / is an H-periodic function, then it follows from (|3.4p that 

6 E cj n) /(i) = 2[/(l,-l,0) + /(l,0,-l) + /(0,l ) -l) + /(-l,l,0)] 

j6H*\H„ 

i<i<«-i 

= 3 £ + + i-l,l)] 

l<j<n-l 

+ 4[/(0,-l,l) + /(-l,0,l)] 

= 6 £ (i- c f 3 )/(i) = 6 E( 1 - c f ) )^- 

jeH„\H*_ 1 jeH„ 
Consequently, we have 

E c S ,l) /(^)-E c S" ) /(i)+ E ^f(i) 

=E c r ) /(i)+E( i -4" ) ) /( ^-5: /( ^ 

Replacing / by fg, we have proved that (f,g) n = (/, g)* n f° r -ff-periodic functions, 
which clearly applies to functions in 7i n . □ 
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3.3. Interpolation on the hexagon. First we note that by ()3.2|) . Theorem 12.61 
when restricted to the hexagon domain becomes the following: 



Proposition 3.4. For n > 0, define 



3n 2 



for f G (7(0). TTien X n / G 7Y„ and 

2n/(^) = /(i) 



Vje 



Furthermore, and perhaps much more interesting, is another interpolation oper- 
ator that works for almost all points in H* . First we need a few more definitions. 
We denote by <9H* := H* \ H° , the set of points in H* that are on the boundary 
of the hexagon {t : — n < ti,t2,h < n}. We further divide effl* as U , where 
consists of the six vertices of dW n and H° consists of the other points in <9H„, 
i.e., the edge points. For j G 9H* , we define 



n n 



— mod 3 



Because of the tiling property of the hexagon, if j G H° then Sj contains two points, 
j and j* G H° , where j* is on the opposite edge, relative to j, of the hexagon; while 
if j G H^, then Sj contains three vertices, j and its rotations by the angles 2ir/3 
and Att/3. 

Theorem 3.5. For n > and f G C(fi), define 



f(i)kn(t) 



where 
(3.7) 



4,«(t)=*;(t-i) and <f:(t)-^E c j"Vj(t) 



jGH* 



Thenl*f€H* n := {0j : j G H;} and 



(3.8) 



J G 



Furthermore, $^(t) is a real function and it has a compact formula 



(3.9) **(t) 



1 

3n 2 



(e n (t) - e n _ 2 (t)) 



(tl - t 2 ) + COS ^(i 2 - i 3 ) + COS 2SZL(i 3 - 



where 
(3.10) 



e„(t) 



Sin Mj^g sin (n+D^2-t3) sin (n+l).(t 3 - tl ) 



sin » 3 ; sin 



3 3 



Proof. First let k G H„. From the definition of the inner product 



it follows 



that if j G H° , then ^j, n (r) = (0k, <f>j)n = ^kj- If j G <9H* , then one of the points in 
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Sj will be in H n , call it j*, then £j, n (^) = ^k,j* ■ As the role of j and k is symmetric 
in £j, n (^), this covers all cases. Putting these cases together, we have proved that 

1 if ^ = i mod 3 



(3.11) =$*J- ~ L ) = < ' " n ~ n 

V ' 3 ' yn> nVn n ' |0 otherwise 



from which the interpolation 

>0 



follows immediately. 



By the definition of cj , it is easy to see that 



where D n is the analogue of the Dirichlet kernel defined by 
(3.12) D B (t) := J2 

By the definitions of 0k and we have 

tfj (t) = 2 [cos 2as: ( t ! _ t 2 ) + cos *f- (t 2 - t 3 ) + cos 2aE (t 3 _ i!)] , 

Jen- 

so that we only have to show that D n = 0„ — O n -i, which in fact has already ap- 
peared in [21j . Since an intermediate step is needed later, we outline the proof of the 
identity below. The hexagon domain can be partitioned into three parallelograms 
as shown in Figure 4. 




Figure 4. The hexagon partitioned into three parallelograms. 
Accordingly, we split the sum over H* as 

— In — In — In 

(3.13) D n (t) = i+ £ £>(t) + E E Z>(*)- 

ia=— nj'i=0 J3=— 7ij2=0 3i=— "33=0 

Using the fact that t\ + t 2 + £3 = and ji + j% + j'3 = 0, each of the three sums can 
be summed up explicitly. For example, 
-1 n 

(3.14) Yl X>(*) 

32=— "31=0 

e _ ¥( „+l) fe _ t2 ) gin _ ^ e _ ¥ „ (tl _ t3 ) gin («+!)£ _ 



sinf (i a - t 3 ) 



sin — t 3 
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The formula D n = 0„ — n _i comes from summing up these terms and carrying 
out simplification using the fact that if t\ + t 2 + £3 = 0, then 

sin 2t 1 + sin 2t 2 + sin 2t^ = — 4 sin t\ sin t 2 sin £3 , 

^ cos2ii + cos2i 2 + cos2i3 = 4cosii cosi 2 cosi3 — 1, 

which can be easily verified. □ 



The compact formula of the interpolation function allows us to estimate the 
operator norm of I* : C(f2) i— > C(f2) in the uniform norm, which is often referred 
to as the Lebesgue constant. 

Theorem 3.6. Let ||2,*||oo denote the operator norm ofT^ : C(Q) 1— » C(fi). Then 
there is a constant c, independent of n, such that 



< c(logn) z 



Proof. A standard procedure shows that 



max 
ten 



E i*;(t-i)i- 



By 



and the fact that D n = <d n — O n _x, it is sufficient to show that 
1 



— — max 

3n 2 ten 



E l^n(t-i)|<c(logn) 2 



n > 0. 



By (|3.13p and (|3.14|) . the estimate can be further reduced to show that 



(3.16) J„(t) := — £ 



sm^f (t 2 - t 3 - &^-) sin 



(n+lW 



jeK* 



and two other similar sums obtained by permuting ti,t%, £3 are bounded by c(logn) 2 
for all t G fi. The three sums are similar, we only work with (|3.16[) . 

Fix a t G 51. If s\ + s 2 + S3 = 0, then the equations s 2 — s 3 = t 2 — t 3 and 
si — s 3 = t\ — £3 has a unique solution s 2 — ^2 and S3 = t 3 . In particular, there 
can be at most one j such that the denominator of F n (t — -) is zero. For t G O, 
setting si = (tx - * 8 )/3 = (2ix + i 2 )/3 and s 2 = (i 2 - i 3 )/3 = (*i + 2t 3 )/3, it 
is easy to see that —1 < si,s 2 < 1. The same consideration also shows that 
— 3n < ji — j3, j 2 — J3 < 3n for j G H* . Consequently, enlarging the set over which 
the summation is taken, it follows that 



1 

3r? 



3n 

E 



- 3 n 



sinn7r(s 2 



k-2 \ 
3nJ 



sin^(s 2 -fj) 



3 71 

E 



fei= 



-3n 



sin(n + 1)tt(si - 



sin7r(si - |M 

V x 3n I 



For s G [— 1, 1], a well-known procedure in one variable shows that 

3n 



IX — 

1.11 n ' 



max 

«e[- 



fc— — 3n 



3n ' 



sin7r(s- 3^) 



< clogn, 



from which the stated results follows immediately. 



□ 
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4. Discrete Fourier analysis on the triangle 

Working with functions invariant under the isometries of the hexagon lattice, we 
can carry out a Fourier analysis on the equilateral triangle based on the analysis 
on the hexagon. 

4.1. Generalized sine and cosine functions. The group Q of isometries of the 
hexagon lattice is generated by the reflections in the edges of the equilateral trian- 
gles inside it, which is the reflection group A 2 . In homogeneous coordinates, the 
three reflections o~\ , a 2 , 03 are defined by 

t<Ti := — (tj., *3, t 2 ), ta 2 := -it 2 , h,t 3 ), ta 3 := -(t 3 , t 2 , £i). 

Because of the relations o~ 3 = a\o 2 <j\ = a 2 G\a 2 , the group Q is given by 

Q = {!,0'l,0'2,0'3,f7iCr2,C r 2a-i}. 

For a function / in homogeneous coordinates, the action of the group Q on / is 
defined by c/(t) = /(to - ), cr G Q. Following [9J, we call a function / invariant 
under Q if erf = / for all a € Q, and call it anti- invariant under Q if of = p(o~)f 
for all a G, where p(o) = —1 if a = o\, a 2l o 3 , and p(c) = 1 for other elements in 
Q. The following proposition is easy to verify ([SJ)- 

Proposition 4.1. Define two operators V + and V~ acting on f(t) by 

(4.1) V ± f(t) = i [/(t) + /(ta 1( r 2 ) + /(ttr 3CTl ) ± /(t CTl ) ± /(ta 2 ) ± /(t<7 3 )] . 

Then the operators T J+ and V are projections from the class of H-periodic func- 
tions onto the class of invariant, respectively anti-invariant functions. 

Recall that ^k(t) = e 3 . Following [21] . we shall call the functions 

TCk(t) := "P + (/>k(t) =- [^An^a.fcaC*) +^fc 2 ,fc3,fei(t) +0k3,fe I ,ka( t ) 
D 

+ ^-fel,-fe 3 ,-fe2( t ) + ^-fc3,-fcl,-fcs(*) + ( /'-fe3,-fc2,-fcl( t )] ) 

TSk(t) := -'P _ 0k(t) =— [0/ Sl ,fc 2 ,fe 3 (t) + 0fc 2 ,fc 3 ,fci (t) + <j>k 3 Mte{*) 
1 01 

a generalized cosine and a generalized sine, respectively, where the second equal 
sign follows from the fact that (pk(ta) = <?W(t) for every <r E Q. It follows that 
TCk and TSk are invariant and anti-invariant, respectively. 

These two functions appear as eigenfunctions of the Laplacian on the triangle 
and they have been studied in [SJ fTS] EH E3 E3 [23 HZj under various coordinate 
systems. In particular, they are shown to share many properties of the classical 
sine and cosine functions in homogeneous coordinates in |2TJ [22] . Some of those 
properties that we shall need will be recalled below. The advantage of homogeneous 
coordinates lies in the symmetry of the formulas. For example, the equations t\ + 
t-2 + t 3 = and k\ + k 2 + k 3 = can be used to write TCk and TSk in several useful 
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forms, such as 

(4.2) TC k (t) 

(4.3) TS k (t) 



l r 

3 



e 3 



(fe2-fc 3 )(t2-*3) 



cos kinti 



+e f(fa-k3)(*3-ti) cosfci7ri 2 + e T(^-fe)(*i-*2) cosknrt 3 



l r 



e 3 



(fca-fc 3 )(t2-ts) 



sin kiirti 



and similar formulas derived from the permutations of t\,t 2 ,t 3 . In particular, it 
follows from (|4.3[) that TSk(t) = whenever k contains one zero component. 

For invariant functions, we can use symmetry to translate the results on the 
hexagon to one of its six equilateral triangles. We shall choose the triangle as 

(4.4) A :={(h,t 2 ,t 3 ) : h + t 2 + t 3 = 0, < t u t 2 , -t 3 < 1} 

={(h,t 2 ):h,t 2 >o, h+t 2 <i}. 

The region A and its relative position in the hexagon is depicted in Figure 5. 




Figure 5. The equilateral triangle A in the hexagon. 

The inner product on A is defined by 

(/,<?> A := r^rjj(t)W)dt = 1 jj{h,t 2 )9{h,t 2 )dtidt 2 . 

If fg is invariant under Q, then it is easy to see that (f,g)H = (/, <?)a- 

When TCk are restricted to the triangle A, we only need to consider a subset of 
k€l. In fact, since TCko-(t) = TCk(tcr) = TCk(t) for t e A and a € Q, we can 
restrict k to the index set 

(4.5) A := {k e H : fci > 0,fc 2 > 0,fc 3 < 0}. 

We define the space T as the collection of finite linear combinations of elements in 
{TC k , TS k : k S A}; that is, 

(4.6) T - I KTC k + fo k TSk) : a kl b k e C, H C A, #S < 00 1 . 
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Since TSk(t) = whenever k has a zero component, TSk(t) arc defined only for 
k e A° , where 

A° := {k e H : fci > 0, k 2 >Q,k 3 < 0}, 
which is the set of the interior points of A. These functions are orthogonal on A. 

Proposition 4.2. For k, j 6 A, 

fl, k = 0, 

(4.7) (TC k ,TCj) A =^.jU, keA\A°,k^0, 

U> keA °> 

and for k, j G A° . 

(4-8) (TS k ,TSj) A = ^ kJ . 

This follows easily from (/,<?)# — (/j<7)a, which holds if fg is invariant, the 
fact that TC k is invariant and TS k is anti-invariant, and the orthogonality of </> k in 
(|3.3p . The norms of TC k are divided into three cases, since TCo(t) = 1 and 

TCk(t) = - [0fc 1 ,fc 2 ,fc 3 (t) + <pk 2 ,k 3 ,ki (t) + ^ks.MM^)} > k 1 k 2 h = 0. 

4.2. Discrete inner product on the triangle. Using the fact that TC k and TS k 
are invariant and anti-invariant under Q and the orthogonality of </> k with respect 
to the symmetric inner product (•,•)£, we can deduce a discrete orthogonality for 
the generalized cosine and sine functions. For this purpose, we define 

A„ := {k e A : -k 3 < n} = {(h,k 2 ) € % 2 ■ h > 0, k 2 > 0,ki + k 2 < n}, 

and denote by A°, A°, A^ the set of points in the interior of A„, on the boundary 
of A„ but not vertices, and at vertices, respectively; that is, 

A° := {k G A„ : ki > 0, k 2 > 0, -fc 3 < n}, 

A v n :={(n,0,-n),(P,n,-n), (0,0,0)}, 

A^ := {k e A„ : k i A° U A v }. 
We then define by TC n and TS n the subspace of T defined by 

TC n = span{TC k : k S A n } and TS n — span{TS k : k S A° }, 
respectively. 

Theorem 4.3. Let the discrete inner product (•, ')A,n be defined by 
where 









je A° 


(4.9) 


A<»> := . 




je A c 








je A v 



Then 



(4.io) (f,g)A = (f,g)A,n, f,g^TC, 
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Proof. We shall deduce the result from the discrete inner product (•, •)* defined in 
(|3.6p . Let ]Q denote the orbit of j under G, which is the set {jo- : a 6 Q}. It is easy 
to see that for j 6 H, 

{6 if hhh ^ 
3 if hhk = but (0,0,0) . 
1 if j = (0,0,0) 

Assume that / is invariant. Using (|4.11[) . we see that 

E ^ n) /(i) = 6E^ n) /e+6 E c l n) /(i) 
= E A j n) /(i), 

and 

E cl n) /(|) = /(o,o,o)+3 4"^) + 3 E c j n) /(|) 
= E A i n) /(i)- 

jeA„ jij 2 j3=o 

Replacing / by /<? and adding these two sums, we have proved that {f,g)* n = 
(/,<?) A, n for invariant fg, which applies to /, g G TC n _i. Thus, the stated result 
follows immediately from Theorem 13.31 and the fact that {f : g)n = (/, s)a- D 

The proof of the above proposition also applies to /, g S 7S n as /g is invariant 
if both / and g are anti-invariant. Noticing also that TSj(-^) = when j e A c , and 
hence we conclude the following result. 

Proposition 4.4. Let £/ie discrete inner product (•, -)A°,n oe defined by 



jeA° 

TTien 

(4.12) (f,g)A = (f,g)A°, n , f,g^rs n . 

The relation (|4.10p can be regarded as a cubature formula 

which is exact for / = g/i, £ TC n . It turns out, however, that much more is 
true. We state the following result in terms of the triangle in (ti,ta)- 

Theorem 4.5. For n > 0, the cubature formula 

(4.13) 2 / ffr^dtrdh = — 2 E E ^/(^f ) 

JA ii=oj 2 =o 

is exact /or aZZ / S TC2n-i- 
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(0,n,-n) 

Q 

<S 

O • (0,n,-n) 

«! : ;»; • o 

• • • * 
• • * • 
• • » • • 

o • » • • • * • » • 

«' V V V V • A ' ?■ ' ' 



O Q V ,V ■/ ' ' ' 

(0,0,0) (n,0,-n) (0,0,0) (n,0,-n) 

Figure 6. Left: A„ in which •, o and o correspond to points of 
different coefficients in (•,•)«■ Right: • indicates points in A°. 

In order to establish such a result, we will need to study the structure of TC„, 
which, and the proof of the theorem, will be given in Section [5j One may note that 
the formulation of the result resembles a Gaussian quadrature. There is indeed a 
connection which will also be explored in Section [5l 

4.3. Interpolation on the triangle. One way to deduce the result on interpo- 
lation is by making use of the orthogonality of generalized trigonometric functions 
with respect to the discrete inner product, as shown in the first theorem below. 
Recall the operator V defined in (|4.ip . 



Theorem 4.6. For n > and f £ C(A) deft 



ne 



C n f(t) := £ /(i)4° n (t), q. n {t) = -j £ TS k (t)TS k (l). 

jeA° ke A £ 
Then C n is the unique function in TS n that satisfies 

£»/(£) = /(£), jeA°. 

Furthermore, the fundamental interpolation function £? n is real and satisfies 

<•?,»(*) = [©«-l(t - |) - ©n-2(t - i 

where means that the operator V~ is acting on the the variable t and 0„ is 
defined in (|3.10[) . 

Proof. By (g^J) and (|4~T2j) . (TSj, TS k ) A °,„ = g^kj, which shows that ^° n (£) = S k<J 
and verifies the interpolation condition. Furthermore, since TS k (J-) = TSj(^) = 
whenever k € A„ \ A°, and TSj(t) = whenever J1J2J3 = 0, it follows from the 
definition of TSj that 

' k £H° 



DISCRETE FOURIER ANALYSIS ON HEXAGON AND TRIANGLE 



19 



where D n is defined in (13.12|) . To complete the proof we use the fact that if F is 
an invariant function, then 

VtV~F(t - s) = 7> t -F(t - a), 

which can be easily verified by using the fact that the elements of the group Q 
satisfy aj = 1, i = 1, 2, 3 and 03 = (T\(t<lO\ = 02 01 02. Q 

We can derive the result of interpolation on the A n using the same approach, 
but it is more illustrative to derive it from the interpolation on the hexagon given 
in Theorem 13.51 

Theorem 4.7. For n > and f G C(A) define 

(4.14) C;j(t) := ]T /(£)#„(t), # n (t) = A<"V + ^„(t), 

where Aj and fj.„ are defined in (|4.9p and (13 . T[) . respectively. Then L* n is the 
unique function in TC n that satisfies 

Furthermore, the fundamental interpolation function £- } _ n is given by 

\v + {@ n {t - i) - 9 n _ 2 (t - I)) - ft {TC„, ,_„(t)TC„,o,-„(i)} • 

Proof. Throughout this proof, write £j = ly n . Recall the equation (|3.11[) . which 
states that ) = 1 if k = j mod 3 and ^j(jjj) = otherwise. Let k £ A. If 
j G A°, then we have F + £j(±) = = [A, (n V^kj- If j G A c and ,hj 2j3 ? 0, 

then V+l } (£) = + <](£)] = |*kj = [Aj^J-Mtj, where fc* is defined 

as in the proof of Theorem 13.51 If j £ A c and j^HH — 0j then 'P + ^j(-) = 
= [^ n) ] _1 ^,j since j G H°. If j G A v and j ^ (0,0,0), then = 
5 Sig5- ^i(jj) = ^j k = [Aj'^l ^^k.jj where <Sj is defined as in the proof of Theorem 
1331 Finally, if j = (0,0,0), V+i^) = 5 k j = [A^'V^kj- Putting these together, 
we have proved that = <5k.j, which verifies the interpolation condition. 

The explicit formula of £^ follows from (|3.9[) and the fact that 

v+ e M*-i)=-p + E mvmS= E TC kW«5 

keH^ keH^ kGH^ 

= 3TC n) 0,- n (t)TC ni o,-n( — ) + 3TCo,n,-ri(t)TCo,n,-n(^), 

as well as the fact that TCo, n ,- n (t) = TC ni o,-n(t). □ 

As a consequence of Theorem 13.61 we immediately have the following: 

Theorem 4.8. Let \\C n \\ and ||£* || denote the operator norm of C n : C(A) 1— * C(A) 
and C* n : C(A) <— > C(A), respectively. Then there is a constant c, independent of 
n, such that 

\\Cn\\ < c(logn) 2 and \\C*J < c(logn) 2 . 
In particular, if f £ C 1 (0), the class of functions having continuous derivatives on 
n, then both C n f and L* n f converge uniformly to f on ft. 
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This theorem shows a sharp contrast between interpolation by trigonometric 
functions on A„ or A° and the interpolation by algebraic polynomials on A n . From 
the result in one variable [17j . the Lebesgue constant for algebraic polynomial in- 
terpolation grows exponentially as n — > oo and, as a consequence, the convergence 
does not hold in C 1 (fi). 

5. Generalized Chebyshev polynomials and their zeros 

Just like in the classical case of one variable, we can use the generalized sine and 
cosine functions to define analogues of Chebyshev polynomials of the first and the 
second kind, respectively, which are orthogonal polynomials of the two variables. 

5.1. Generalized Chebyshev polynomials. To see the polynomial structure 
among the generalized trigonometric functions, we make a change of variables. 
Denote 

(5.1) z := TCo,i,-i(t) = !0o,i,-i(t) + 0i,_i,o(t) + <^-i,o,i(t)], 

whose real part and imaginary part, after simplification using (|3.15p . are 

x = JRz — | cos j(t 2 — ti) cos ^(^3 — t 2 ) cos j(h — t 3 ) — i, 

y = 3z = | sin|(t 2 - ii)sin|(t 3 - t 2 )sin|(t 1 - t 3 ), 

respectively. If we change variables (tijtg) l— * ( x >y)> then the region A is mapped 
onto the region A* bounded by Steiner's hypocycloid, 

(5.3) A* = {{x, y) : ~3{x 2 + y 2 + l) 2 + 8{x 3 - 3xy 2 ) + 4 > 0}. 



Figure 7. The region A* bounded by Steiner's hypocycloid. 

Definition 5.1. Under the change of variables ()5. 1|) . define 
T™(z, z) : = TC fe , m _ fc ,_ m (t), < k < m, 

UJ?(z, z) : = TS fc+i^-fc+i,-rn- 2 (t) ^ o < fc < m. 
TSi,i,_2(t) 

We will also use the notation T™(t) and U™(t) when it is more convenient to work 
with t variable. 
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These functions are in fact polynomials of degree m in the z and z variables since 
they both satisfy a simple recursive relation. 

Proposition 5.2. Let Pig 1 denote either T£ l or . Then 

(5.4) 



P%_ k (z,z)=P?(z,z), 0<k<m, 
and they satisfy the recursion relation 

(5.5) Pl n+1 (z, z) = 3zP™(z, z) - P™ +1 (z, z) - P£i l (z, z) 

for < k < m and m > 1, where we use 



T^(z,z)^Tr +1 (z,z), T™ +1 (z,z) 



rpm-\-l 



U? 1 (z,z) = 0, U%~ 1 (z,z) = 0. 
In particular, we have 

T a( z , z) = !> T o( z , z ) = z -, T i( z > z ) = f > 



U°(z,z) 



Uq(z,z)=3z, U{(z,z) = 3z. 



The recursion relation (|5.5j) follows from a straightforward computation. To- 
gether, (|5.5[) and (|5 .4|) determine all P™ recursively, which show that both T™ and 
U™ are polynomials of degree m in z and z. Furthermore, each family inherits 
an orthogonal relation from the generalized trigonometric functions, so that they 
are orthogonal polynomials of two variables. They were studied by Koornwindcr 
([SIHO])) who called them the generalized Chebyshev polynomials of the first and 
the second kind, respectively. Moreover, they are special cases of orthogonal poly- 
nomials with respect to the weighted inner product 



(f,g)w a = c a / f{x,y)g{x,y)w a (x,y)dxdy 1 



where w a is the weight function defined in terms of the Jacobian of the change of 
variables (|5.1[) . 

2a 

1 1\ :■ in lit 

w a (x,y) 



d(x,y) 


2a r 


d(h,t 2 ) 





16 2 . . . . 

— 7r sm 7Tix sm 7ri2 sin7r(ti + t 2 ) 
4 Q 



= [~3(x* + if + If + 8(x 3 - 3^ 2 ) + 4] a 

and c a is a normalization constant, c a :— 1/ J A w a (x,y)dxdy. We have, in partic- 
ular, that c_i = 2 and ci = 81/(87r 4 ). Since the change of variables (|5.ip shows 
immediately that 

T^l J f(t)dt = c_i J f(x,y)w_i(x,y)dxdy, 

it follows from the orthogonality of TCj and TSj that T™ and U™ are orthog- 
onal polynomials with respect to w_i and wi, respectively. Furthermore, from 
Proposition 14. 21 we have 



(5.6) (T™ ,Tp) w _ 1 =d^5k,j, 



df ■= 



m = k = 0, 
k = or k ■ 
1 < k < m 



m, m > 0, 
l,m > 0, 
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and 

(5.7) (U?, U™) W1 = U kJ , < j, k < m. 

2 

Let T™ = ^[S^T™ and {7™ = ^/6U™. These are orthonormal polynomials. We 
introduce the notation 

(5.8) F m = {P ™,P?,...,P%} 

for either P™ = TJJ 1 or UJJ 1 . The polynomials in P„ consist of an orthonormal basis 
for the space V n of orthogonal polynomials of degree m. Evidently both T™ and 
U™ are complex polynomials. As a consequence of f|5 .4[) . however, it is easy to 
derive real orthogonal polynomials with respect to w ± i(x,y) on A*. In fact, let 
P m be as in (|5.8[> . Then the polynomials in the set P m defined by 

vL-i(^v) :={^«{i , * ,n - 1 (^^)},V^>{Jf m - 1 (^2)},0<fe< Lm-lJ} 

PL(z,y) :={^ m , V^P^,*)}, V23{P|'"(z,z)},0 < k < |m|} 

are real orthonormal polynomials and form a basis for V n with n = 2m — 1 or 2m, 
respectively. 

5.2. Zeros of Chebyshev polynomials and Gaussian cubature formula. Let 

w be a nonnegative weight function defined on a compact set fl in ]R 2 . A cubature 
formula of degree 2n — 1 is a sum of point evaluations that satisfies 

A 

/ f(x,y)w{x,y)dxdy = V \ j f{x j ,y j ), Xj G M 

for every / G n| n _ 1 , where IT 2 denote the space of polynomials of degree at most 
n in two variables. The points (xj,yj) are called nodes of the cubature and the 
numbers A; are called weights. It is well-known that a cubature formula of degree 
2n — 1 exists only if N > dimll 2 ^. A cubature that attains such a lower bound is 
called a Gaussian cubature. 

Let P n := {PJ} : < k < n} denote a basis of orthonormal polynomials of degree 
n with respect to w(x, y)dxdy. It is known that a Gaussian cubature exists if and 
only if P„ (every element in P„) has diml! 2 ^ many common zeros, which are 
necessarily real and distinct. Furthermore, the weight Xj in a Gaussian cubature is 
given by 

n k 

(5.10) X 3 = [K^x^yj)]- 1 , K n (x,y) := £ £ Pf(x)Pf(y). 

k=0 j=0 

The definition of K n (x,y) shows that it is the reproducing kernel of IT 2 , in L 2 (fi), 
so that it is independent of the choice of the basis. 

Unlike the case of one variable, however, Gaussian cubatures do not exist in 
general. We refer to [H HH [5D] for results and further discussions. It is sufficient 
to say that Gaussian cubatures are rare. In fact, the family of weight functions 
identified in [18] for which the Gaussian cubature exist for all n remains the only 
known example up to now. 

We now test the above theory on the generalized Chebyshev polynomials. Be- 
cause of the relation ()5.9|1 , we can work with the zeros of the complex polynomials 
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in P n . Furthermore, it follows from (15. 4| that the reproducing kernel 

n n k 

K n (z, w) := Ki^ z)¥ k (w, = E E P /( z ' ^ P H W > 
fc=o fe=o j=o 

is real and, by (|5.9p , is indeed the reproducing kernel of IT^ . 

For the Chebyshev polynomials of the first kind, the polynomials in P n do not 
have common zeros in general. For example, in the case of m = 2, we have 
Tq(z, z) = 3z 2 — 2z and T 2 (z,z) = (3zz — l)/2, so that the three real orthogo- 
nal polynomials of degree 2 can be given by 

{3(x 2 - y 2 ) - 2x, 6xy + 2y, 3{x 2 + y 2 ) - l} , 

which has no common zero at all. Therefore, there is no Gaussian cubature for the 
weigh function w_i on A*. 

The Chebyshev polynomials of the second kind, on the other hand, does have the 
maximum number of common zeros, so that a Gaussian cubature exists for wi(t). 
Let 

X n := A° +2 = {(^s, > 1, k 2 > 1, h + k 2 < n + 1}, 

which is in the interior of A, and let X* denote the image of X n in A* under the 
mapping (|5.2|) . 

Theorem 5.3. For the weight function wi on A* a Gaussian cubature formula 
exists; that is, 



(5.11) ci / f{x,y)w k {x,y)dxdy=Y J E ^/(^^ V/eH 

JA * Ji=l 32 = 1 

where j = (ji, j 2 ), 



2n-l) 



Proof. The numerator of UJ}(t) is TSfc+i.„_fc+i._ n _2(t), which is zero at all 
according to (|4.3[) . The first equation of (|3.15[) shows that the denominator of 

um is 

4 

(5.12) TSi j i._2(t) = — sin7r^i sin7ri2 sin7r^3, 

3 

which vanishes at the integer points on the boundary point of A. Consequently, 
Uj}(t) does not vanish on the boundary of A but vanishes on X n . It is easy to see 
that the cardinality of X n is dim 11^ Consequently, the corresponding P n has 
dimil 2 ^ common zeros in X*, so that the Gaussian cubature with respect to wi 

exists, which takes the form (|5.11[) . To obtain the formula for fJ^, we note that 
the change of variables (|5.2[) shows that (|5.1ip is the same as 

± J a f(t) [TSr^it)} 2 dt = £ [ TS M.- 2 (^)] 2 /(^), 

where the last step follows from the fact that TSi.i^Ct) vanishes on the boundary 
of A. This is exactly the cubature (|4.13[) applied to the function /(t) [TSi.i i _2(t)] 2 . 

□ 
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We note that wi is only the second example of a weight function for which a 
Gaussian cubature exists for all n. Furthermore, the example of w_i shows that the 
existence of Gaussian cubature does not hold for all weight functions in the family 
of w a . In 18J, a family of weight functions depending on three parameters that 
are on a region bounded by a parabola and two lines is shown to ensure Gaussian 
cubature for all permissible parameters. 

The set of nodes of a Gaussian cubature is poised for polynomial interpolation, 
which shows that there is a unique polynomial P in H^ l _ 1 such that = /(£) 
for all £ € X*, where / is an arbitrary function on X*. In fact the general theory 
(cf. [4]) immediately gives the following result: 

Proposition 5.4. The unique interpolation polynomial of degree n on X* is given 
by 

n+l n+1— ji 

(5-13) £*/(*)=£ E /(M°*n(t,£), 

jl=l J2=l 

where x = (x\, x%) and z(t) = x% + ix 2 , and C n f is of total degree n — 1 in x. 

Under the mapping (|5.2[) , the above interpolation corresponds to a trigonometric 
interpolation on A° +2 . Since = ( n + 2 ) 2 |~TSi ; i i _2(^) | 2 , it follows readily that 

Kn{t > " } - fr + W TS ^-*(n) TS M ,_ 2 (t) ' 

where j = (ji, j2, — ji —32)- Consequently, the interpolation (|5.13[) is closely related 
to the trigonometric interpolation on the triangle given in Theorem 14.61 



5.3. Cubature formula and Chebyshev polynomials of the first kind. As 

discussed in the previous subsection, Gaussian cubature does not exist for w_i on 
A*. It turns out, however, that there is a Gauss-Lobatto type cubature for this 
weight function whose nodes are the points in the set 

Y n := A„ = ^)\k t >0,k 2 > 0, h + k 2 < n}, 

which is located inside A. Again we let Y* denote the image of Y n in A* under the 
mapping (|5.2|) . 

The cardinality of Y n is dim 11^. According to a characterization of cubature 
formula in the language of ideals and varieties ( 25 ), there will be a cubature 
formula of degree 2n — 1 based on Y* if there are n + 2 linearly independent 
polynomials of degree n+l that vanish on Y*. In other words, Y* is the variety 
of an ideal generated by n + 2 algebraic polynomials of degree n+l, which are 
necessarily quasi-orthogonal in the sense that they are orthogonal to all polynomials 
of degree n — 2. Such an ideal and a basis are described in the following proposition. 

Proposition 5.5. The set Y* is the variety of the polynomial ideal 

(5.14) (T " +1 -Tr, T^-T^ (l<fc<n), T^-T^). 
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Proof. The mapping (|5.2| allows us to work with TJ}(t) on Y n . Using the explicit 
formula for TCfc n-fc.-n in (14.21) , we have 



1 r hi 

e 3 



(n-2fc)(t 2 -t 3 ) 



cos Tint -i 



+e ^(n-2fc)(t,- tl ) cogn7rt2 + ^ (n-^^-t,) CQSmrt3 

Since n + 1 — 2k = n — 1 — 2(fc — 1), the above formula implies that 



T fe n+1 (t)-T^(t) = ^ 



e ^(n-2fc)(t 2 -t 3 ) sm7r ^ sinm7r ^ 
+ (n-2*)(*3-ti) gin7r ^ sinn7rt2 



-e 3 



Ht(„_2fe)(ti-t 2 ) 



sin 7ri3 sin mrtz 



from which it follows immediately that T^ +1 — T2\ vanishes on Y n for 1 < fc < n. 
Furthermore, a tedious computation shows that 



nfji 32. ~Ji~32 \ _ 



fc V n ' n 



-(n+l)ji+j 2 



(-2 + e 2inj2 + e 2i7r(j ' 1+J ' 2) ) 



+e 2 ^f (_ 2e 2 ^2 + e 2i7r(j 1+ j 2 ) + !) + e 2l,r2i ^ 2a (-2e 2i ^ 1 + e 2i7rj2 + 1) 



which is evidently zero whenever j\ and j% are integers. Finally, we note that 
T^Xl - Tn-\ is the conjugate of T ,l+1 - T x n . This completes the proof. □ 

As mentioned before, this proposition already implies that a cubature formula 
of degree 2n — 1 based on the nodes of Y* exists. 

Theorem 5.6. For the weight function w_i on A* the cubature formula 
(5.15) 

c_i / f(x, y)w_i(x, y)dxdy = £ £ £), V/ € n an _!, 

A * U fc 1= 0fc 2 =0 



holds, where X% = A^_ fc _ n im£/i Aj™" 1 given by (|4.9p . 

In fact, the change of variables back to A shows that the cubature (|5 . 1 5[) is 
exactly of the form (|4.13|) , which gives the values of X^ . In particular, this gives 
a proof of Theorem 14.51 

Another way of determining the cubature weight Ajj, is to work with the repro- 
ducing kernel and the polynomial ideal. This approach also gives another way of 
determining the formula for polynomial interpolation. Furthermore, it fits into a 
general framework (see |24) and the references therein), and the intermediate results 
are of independent interest. In the rest of this section, we follow through with this 
approach. 

The recursive relation (|5.5|) for T™ can be conveniently written in matrix form. 
Recall that for T™, the orthonormal polynomials are 



An) 



P™ = V /3T "\ Pf = x/6Tf (l<fc<m-l), P% = V3T% 
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for m > and Pg = 1. In the following we treat P m defined in (|5.8j) as a column 
vector. The relation (|5.4| shows an important relation 

— \ • 1 

(5.16) P m = J m -j_]P m , J m = 

.i OJ 

Using these notations, we have the following three-term relations: 
Proposition 5.7. Define P_i = 0. For m > 0, 

— ^mlPm+1 ~l~ -^m^m -\- O m P m _i, 
zP m = Cra+l^m+l + ^m^m + "^m— l^m-l i 

where A m , B m and C m are given by 

"0 V2 ... 



(5.17) 



1 



.0 



O o- 



y/2 





o O 





1 



O 



1 

y/2 



and C m — Jm-AJJ, 



\Jm— 1- 

In fact, the first relation is exactly (|5.5[) . To verify the second relation, take the 
complex conjugate of the first equation and then use (|5.16p . 

In the following we slightly abuse the notation by writing P„(t) in place of P„(z), 
when z = z(t) is given by (|5.1[) . Similarly, we write i^T n (t,s) in place of K n (z,w) 
when z = z(t) and w = w(s). Just like the case of real orthogonal polynomials, the 
three-term relation implies a Christoph-Darboux type formula (|5j) 

Proposition 5.8. For m > 0, we have 

(5.18) (z - w)K n (z,w) := P^ +1 (z)< r P^J - W*(z)C% +1 F n+1 (w). 

Next we incorporate the information on the ideal in (I5.14|) into the kernel func- 
tion. The n + 2 polynomials in the ideal can be written as a set 

(5.19) Q n+1 = p„+i - r p„ - rxPn-i, 

which we also treat as a column vector below, where 





"0 


7S 


... 0" 


r = 


O 


... o 


... o 







... 





and Ti 





V2 



o 















o 



V2 




Proposition 5.9. For n > 0, define 

$„(t,s) := i [K n (t,s)+K n ^(t,s)} - ~ [r "(t)Ty(i) + 2^t)T n (s)] , 
and we again write <& n (z, w) in place of $„(t, s) for z = z(t) and w = w(s). Then 
(5.20) (z - w)$ n (z, w) = Q« +1 (z)A«M n ¥jw) - P£( 2 )M n C^iQ^H 
where M n = diag{— 1/3, 1/2, . . . , 1/2, — 1/3} is a diagonal matrix. 
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Proof. Inserting P„+i = Q n+1 
(z - w)K n (z,w) 



r P„ + TiP„_i into (PS)) gives 



£ + l(*R r Pn(«0 - P^)C+iQ»+iH 



+ P^rg^ p„h - p t „ r (z)C + ir P„(^) 



Let D n be a n x n matrix. Then similarly we have 



(z - w)F«(z)D n ¥Jw) = Q« +1 {z)A«D n ¥Jw) - W%(z)D n C% +1 Q n+1 {w) 

+ F«(z)(A n r + B n ) tr D n ¥Jw)-F«(z)D n (C% +1 r + B%)¥Jwj 

+ F t ; i _ 1 (z)(A n T 1 + C n f r D n ¥jw) - F t ;(z)D n {C^ +1 T 1 + A^F^Jw). 

In order to prove the stated result, we take the difference of the two equations and 
choose D„ such that 



(5.21) 



(^4„r + B n ) tr D n - D n (c% +1 r + b„) — TqA^ - c* r +1 ri, 

(A n Ti + C n fD n = r*{An, DniC^T, + Al_ x ) = C« +l T x . 



It turns out that the unique solution of the above three equations is given by 
D n = diag{2/3, 1/2, 1/2, . . . , 1/2, 2/3}, which is completely determined by the last 
two equations and then verified to satisfy the first equation. As a consequence of 
this choice of we obtain 



K n (z,w)-F*(z)D n F(w) 



■■(T< +1 (z)A«(I-D n )F n (w) 
F%(z)(I - D n )C% +1 ® n+1 ( W ) 



where I denotes the identity matrix. Setting M n = I — D n we see that the right 
hand side agrees with the right hand side of (|5.20| . Finally we note that 

1, 



Fl(z)D n F n (w) = -Fl{z)F n (w) + [P£(z)P£(w) + P£(z)P<l(w) 

= o l K n(z,w) - K n -i(z,w)] + i 
which completes the proof. 



T£{z)TS{w)+TZ{z)TS{z) 



As a consequence of the above proposition, it is readily seen that ^ n { L ,-) = ( 
for j 7^ k. Furthermore, we note that $„(t,s) is in fact a real valued function 
Consequently, we have proved the following. 

Proposition 5.10. The unique interpolation polynomial on Y* is given by 



n n—]x 



where z(t) = x + iy, and C* n f is of total degree n in x and y. 



In fact, upon changing variables (|5.1[) . this interpolation polynomial is exactly 
the trigonometric interpolation function £*/(t) at (|4. 14|) . Here the uniqueness of 
the interpolation follows from the general theory. Moreover, integrating £*/ with 
respect to w_i over A* yields the cubature formula (|5.15|) with 



(5.22) 



A 



O) 
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which can be used to determine X" if needed. The function K n (t, s), thus < I ) n (t, s), 
enjoys a compact formula, which facilitates the computation of <& n (^,^). The 
compact formula follows from the following identity. 

Proposition 5.11. Let D n be defined in (I3.12|) . Then 

n k 

K n (s,t) = ££fj=(t)lJ( S ) = V+D n (t s). 

fc=0 3=0 

The identity can be established from the definition of and D n as in the proof 
of Theorem l4.6l upon using the fact that V^V+F(t-s) = V+F(t-s). We omit the 
details. From this identity and (|5.22p , we immediately have a compact formula for 
the interpolation polynomial, whose trigonometric counterpart under the change of 
variables (|5.2[) is exactly the function £^ n in Theorem 14.71 
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